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power  density,  dynamically  similar  units  can  be  produced  using  these  laws. 

The  laws  are  developed  following  the  modern  procedure  of  ordering.  Numerical 
examples  are  provided  to  Illustrate  the  procedure.. 
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SECTION  I 
INTRODUCTION 

Designing  MHD  generators  for  high  power  density  Is  currently  limited 
to  a  small  circle  of  practitioners,  away  from  the  main  stream  of  commer¬ 
cial  MHD  activities;  the  planning  Is  usually  done  In  an  ad  hoc  manner, 
for  Individual  cases,  without  stating  clearly  the  approach  followed  or 
the  philosophy  adopted.  To  render  such  planning  more  a  science  than  an 
art.  It  Is  the  purpose  of  this  study  to  formalize  the  steps  to  be  taken 
In  the  preliminary  design  of  a  high  power  density,  open  cycle  segmented 
Faraday  or  diagonal  combustion-driven  MHD  generator.  The  recommended 
steps  are  gathered  here  for  one  purpose:  the  realization  of  maximum 
possible  load  power  per  unit  channel  volume.  Detailed  design  computations 
that  follow  these  preliminary  steps  will  undoubtedly  Indicate  trade-off 
points  for  a  variety  of  meritorious  features. 

In  obtaining  material  from  a  widely  scattered  literature  for  the 
expressed  purpose,  the  author  can  hardly  claim  any  originality.  Although 
critical  comments,  personal  judgments,  minor  discoveries,  slight  exten¬ 
sions  and  small  variations  abound,  this  report  Is  meant  to  be  a  designer's 
guide. 

The  approach  toward  achieving  high  power  density  will  be  Iterative. 
Choice  of  a  magnetic  field  strength  profile  is  regarded  as  the  first 
among  a  series  of  design  decisions  that  will  lead  to  a  practical  system. 
The  magnetic  field  strength  profile  Is  to  be  chosen  from  experience  and 
research  of  the  literature  as  well  as  from  the  following  considerations. 
Use  the  maximum  magnetic  flux  through  the  MHD  flow  consistent  with  all 
the  usual  related  pragmatic  design  factors.  These  magnet  design  factors 
Include  Hall  voltage  breakdown  limit  between  MHD  channel  electrode 
segments,  weight/volume  as  a  function  of  magnetic  field  strength,  and 
cost.  After  going  through  the  preliminary  design,  the  system  must  be 
examined  to  see  If  the  design  Is  practicable.  There  seems  to  be  no 
short-cut  for  arriving  at  the  final  practicable  design.  If  the  magnetic 
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field  profile  Is  not  suitably  defined  In  the  beginning,  the  entire  pro¬ 
cess  must  be  repeated  with  an  adjusted  profile  until  the  system  design 
Is  satisfactory.  (Although  no  attempt  has  been  made  to  computerize  the 
overall  design  procedure,  with  much  effort  one  conceivably  can  build  a 
complex  computerized  design  process  that  starts  with  system  constraints 
and  arrives  at  a  finely  optimized  design.) 

The  design  thus  starts  with  a  prescribed  magnetic  field  strength 
profile  along  the  active  channel  length.  In  Section  II,  the  recommended 
design  process  establishes  the  combustion  chamber  pressure,  after  a  brief 
description  of  optimization  calculations  to  be  done  on  the  inlet  Mach 
number,  seeding  ratio,  and  0/F  (oxygen  to  fuel  ratio).  The  optimization 
Is  based  on  a  semi -empirical  expression  of  the  effective  power  output 
developed  by  Smith  and  Nichols  (Reference  1).  This  new  approach  provides 
a  rational  basis  to  the  design  procedure. 

As  a  result  of  Section  II,  estimate  of  the  realizable  power  density 
emerges,  which  will  yield  the  magnitude  of  the  transverse  dimension  for 
a  fixed  length/diameter  ratio,  for  the  desired  power  output. 

In  Section  III,  It  Is  recommended  that  the  cross-sectional  area 
variation  and  the  length  of  the  generator  be  estimated  on  the  basis  of 
an  isothermal  core  flow.  The  rationale  here  Is  as  follows:  the  electri¬ 
cal  conductivity  of  the  plasma  varies  exponentially  with  temperature; 
therefore,  keeping  the  entire  duct  at  a  uniformly  high  temperature  level 
would  promote  high  power  density.  In  contrast,  a  constant  velocity  de¬ 
sign  would  Incur  heavy  temperature  drop  and  pressure  loss;  the  latter 
would  also  burden  the  diffuser  heavily.  This  preliminary  shape  of  the 
duct  will  eventually  serve  as  the  base  on  which  variations  will  be  made 
on  a  computer  to  accommodate  the  boundary  effects.  The  final  (computer- 
aided)  design,  of  course,  will  not  quite  come  out  Isothermal. 
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Section  III  contains  formulas  developed  for  the  Isothermal  core  flow 
through  a  diagonal  conducting-wall  generator.  To  the  author's  knowledge, 
these  extended  formulas  for  the  case  of  diagonal  conducting  walls  are  new. 

In  Section  IV,  the  preliminary  design  Is  carried  out  along  a 
different  route.  Starting  with  any  well  designed  generator  which  Is 
demonstrably  high  In  power  density,  either  already  In  operation  or  In 
an  advanced  stage  of  planning,  scaling  laws  c?n  be  applied  to  produce 
dynamically  similar  units  for  a  different  power  output,  and/or  a  differ¬ 
ent  magnet  strength,  and/or  a  different  fuel,  etc.  The  procedure  can, 
for  example,  be  used  to  yield  a  small  dynamically  similar  pilot  unit 
which  can  be  tested  before  embarking  on  a  larger-sized  endeavor.  It 
must  be  emphasized  here  that,  If  a  unit  has  a  power  density  which  Is 
maximum  under  the  given  restraints.  Its  dynamically  similar  models  will 
deliver  smaller  power  per  unit  volume  (being  still  "high",  possibly). 
However,  dynamically  similar  models  will  work  with  equal  efficiencies. 

The  key  modeling  parameter  Involved  In  the  scaling  laws  are 
established  In  Section  IV  following  a  modern  procedure  known  as  the 
ordering  process  (see,  e.g..  Chapter  5  of  Reference  2).  Modeling  (or 
scaling)  Is  then  carried  out  In  the  classical  manner  (see,  e.g..  Chapter 
4  of  Reference  3).  Although  the  resulting  laws  are  Identical  with  those 
quoted  In  the  literature  (e.g..  Chapter  7  of  Reference  4),  the  derivation 
presented  In  this  report  may  be  more  convincing. 

Finally,  In  Section  V,  minor  losses  near  the  walls  are  discussed. 

The  discussion  is  brief  since  these  losses  will  also  be  accounted  for 
In  the  final  design— the  computer-aided  simulation  and  selection. 

It  is  hoped  that  designers  will  find  the  recommended  procedure 
helpful  In  providing  Initial  inputs  to  the  design  of  the  magnet,  which 
usually  has  to  be  started  simultaneously  with  that  of  the  generator,  as 
well  as  being  helpful  in  starting  the  sophisticated  computations. 
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SECTION  II 

SELECTION  OF  PRESSURE  LEVEL 

A  modern  and  operational  definition  of  design  Is  "optimization  under 
partially  uncertain  constraints."  With  this  definition  In  mind,  one  may 
state  a  general  design  philosophy  or  approach  In  the  form  of  five  steps. 
Assuming  that  the  operation  Is  not  too  sensitive  to  changes  of  various 
parameters  In  a  rather  large  neighborhood  of  the  optimal  condition,  (1) 
Establish  a  rational  guideline  for  the  optimization  with  respect  to  each 
parameter.  (2)  Optimize  the  object  quantity  with  respect  to  the  para¬ 
meters  one  after  the  other.  (3)  Display  a  number  of  optimal  calculations 
over  a  range  of  uncertain  values  of  the  constraining  parameters.  Select 
a  few  parameter  sets,  exercising  designer's  judgement.  (4)  Trade  off 
(1 .e. ,  deviate  from  the  optimal)  for  other  desired  or  required  character¬ 
istics.  (5)  Model  the  few  cases  which  exhibit  preferred  overall 
characteristics  on  a  computer,  and  make  a  final  optimization  through 
judicious  Iterations.  The  object  quantity  to  be  maximized  In  the  above 
steps  Is  the  load  power  output  per  unit  channel  volume,  w.  At  the  outset, 
the  fuel  used,  the  oxydlzer,  and  the  seeding  material  are  chosen;  but 
these  choices  are  outside  the  scope  of  this  report.  The  seeding  ratio, 

0/F  ratio.  Inlet  Mach  number  Mj,  and  especially  the  combustion  chamber 
pressure  p°  are  among  parameters  to  be  determined.  All  products  of 
combustion  In  this  report  refer  to  hydrocarbon  fuels. 


Effective  power  density  Is  defined  as  the  delivered  load  power 
divided  by  the  active  channel  volume.  As  a  measure  of  the  effective 
power  density  when  loaded  for  maximum  power  under  fixed  propellant  rate, 
the  following  semi-empirical  formulas  from  Smith  and  Nichols  (Reference 
1)  apply: 

Segmented  Faraday— 


Vf ' 


(X/%)o*u*2B*2,  &*  $  1 


e*  >  i 
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Diagonal  Conducting  Wall  — 

■>*  i  \2 


w  .3 
eff 


-1J1.  (V.)aV2B»2,  , 

8*0  +  a)  {1  +  8  2) 


~  ^ - (VOoW2,  3  ,  1 


(1  +  a*2)  (1  +  8*2) 


where  the  asterisk  Is  used  to  Indicate  evaluation  of  quantities  at  a 
certain  reference  point  (the  MHD  channel  Inlet,  for  example),  and  where 

8  =  Hall  parameter 
a  -  conductivity 
u  =  flow  velocity 
B  =  magnetic  field  strength 

a  =  (Hall  field  strength)/ (Faraday  field  strength) 


In  an  empirical  and  approximate  manner,  the  formulas  account  for  the 
Internal  current  leakage  and  electrode  voltage  drops. 


For  a  given  fuel  mixture,  temperature  level,  and  B*.  the  weff  vs. 

8*  curve  (remembering  that  because  of  the  pressure  variation) 

shows  a  trend  as  plotted  In  Figure  1  (the  curve  marked  diagonal  conducting 
walls  being  roughly  for  a=-8,  a  power-maximizing  value).  An  Important 
point  Is  that  the  Illustrated  variation  of  8*  Is  to  be  obtained  solely 
through  varying  the  dimensionless  pressure  p.  (Thus,  at  the  origin, 

8*=  o,  the  pressure  p  Is  infinite).  It  Is  thus  observed  that  wg^f  Is 
maximum  when  8*=1. 

A  rational  guideline  for^the  selection  of  the  pressure  level  now 
emerges:  choose  p  such  that  8=1.  for  the  given  T*  and  B*;  If  practical 
constraints  force  a  deviation,  It  Is  slightly  better  to  make  8*  >  1 
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(rather  than  <1).  The  latter  preference  for  g*>l  Is  based  on  the  obser¬ 
vation  that  the  limiting  Faraday  case  has  a  horizontal  tangent  at  8*  ■  1 
(on  the  6*>1  side)  In  Figure  1. 

In  following  this  guideline,  there  Is  still  the  question  of  where  to 
enforce  It.  Should  the  asterisk  reference  point  refer  to  the  nozzle  exit, 
or  to  where  the  magnetic  field  peaks,  or  should  one  use  some  kind  of 
average  state?  Anticipating  an  Isothermal  (or,  approximately  Isothermal) 
design,  a  will  stay  approximately  the  same  everywhere  In  the  core  flow. 
Also,  from  existing  designs  where  temperature  drop  Is  relatively  slight, 
we  observe  that  the  flow  velocity  hardly  changes  upstream  of  the  peak  of 
the  magnetic  field.  Therefore,  It  seems  desirable  to  select  the  pressure 
value  so  as  to  make  8=1  (or  somewhat  above)  at  the  peak  of  the  magnetic 
field  strength.  Thus,  with  a  computer  print-out  of  8(T,p,B)  In  hand, 
one  will  select  this  pressure  level  accordingly. 

Once  the  pressure  at  the  peak  magnetic  field  strength  position  Is 
chosen,  one  can  estimate  an  additional  percentage  (e.g.,10%)  to  obtain 
p.j  at  the  Inlet.  The  combustion  pressure  can  then  be  calculated  for  a 
given  entrance  Mach  number  Mj,  knowing  the  ratio  y  of  the  specific  heat 
capacities  of  the  plasma.  also  links  the  generator  temperature  with 
the  combustion  gas  temperature  T®. 


As  an  example,  toluene  +  Og  +  Cs2C03  at  the  stoichiometric  0/F  ratio 

Is  treated  In  the  manner  just  described  for  a  peak  magnetic  field  of  4  T. 
The  result  Is  summarized  In  Table  1.  Two  effective  combustion  gas  tempera¬ 
tures  are  employed  In  the  table  since  the  temperature  level  at  the  nozzle 
entrance  may  be  raised  somewhat  by  an  Increase  In  combustion  pressure, 
and/or  an  Improvement  In  the  combustion  chamber  design,  etc. 

It  should  be  emphasized  once  again  that  the  optimal  pressure  level 
Is  chosen  Independently  of  optimization  with  respect  to  the  seeding  ratio, 
0/F  ratio  and  Inlet  Mach  number.  If  all  parameters  are  optimal,  the 
effective  power  density  will  be  the  maximum  of  maxima;  otherwise.  It  will 
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only  be  the  maximum  for  a  given  (non-optlmal )  set  of  parameter  values. 

In  Table  1,  although  the  Mach  number  Is  In  the  optimal  range,  the  seeding 
and  the  0/F  ratios  are  not  necessarily  optimal. 

Alternatively,  with  the  aid  of  a  computer,  the  effective  power  den¬ 
sity  can  be  optimized  with  respect  to  each  of  the  four  parameters,  seed¬ 
ing  ratio,  0/F  ratio.  Inlet  Mach  number,  and  pressure  level.  Essentially, 
the  computer  Is  to  print  out  a  chart  showing  the  variation  of  wg^f  with 
respect  to  these  parameters;  the  optimal  values  are  then  Identified. 

Less  elaborate  study  (Reference  5)  has  shown  that  M^  should  always  be 
around  2  for  products  of  combustion  to  realize  maximum  power  density 
(which  fact  Is  also  borne  out  by  a  detailed  numerical  example  In  (Refer¬ 
ence  1).  Thus,  for  a  preliminary  design,  anticipating  certain  practical 
ranges  of  the  seeding  and  0/F  ratios  (which  should  be,  but  might  not  be, 
around  the  optimal  value)  and  M^  (which  must  be  around  2),  the  pressure 
level  may  as  well  be  selected  by  following  the  suggested  guideline  (i.e., 
to  make  e=l  where  the  magnetic  field  peaks,  to  estimate  p^  by  adding  a 
percentage,  and  to  calculate  p®  based  on  M^)  as  exemplified  In  Table  1. 

Finally,  each  selection  In  Table  1  has  an  anticipated  power  density 
associated  with  It;  we  will  quote  four  numbers  here  as  Illustrations: 

10*  seeding,  M^  *  2.1  — 

T®  «  3100  K:  w  ^60  MW/m3 

T°  -  3400  K:  w  ^320  MW/m3 

30*  seeding,  Mj  ■  2.1  — 

T®  =  3100  K:  w  M25  MW/m3 
T®  «  3400  K:  w  ^620  MW/m3 

From  these  values  the  generator  volume  can  be  estimated  for  a  desired 
power  output.  If  an  empirical  length- to-dlameter  ratio,  e.g.  10,  Is 
adopted,  based  on  a  compromise  between  end  and  wall  effects,  the  dimen¬ 
sions  of  the  generator  can  then  be  estimated. 
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TABLE  1 


OPTIMAL  COMBUSTION-CHAMBER  PRESSURE  FOR  STOI¬ 
CHIOMETRIC  COMBUSTION  OF  TOLUENE  AND  OXYGEN 
SEEDED  WITH  CESIUM  CARBONATE 
(4  TESLA) 


10£  Seeding  By  Wt 

light  Of  Fuel 

30^  Seeding  By  Weight  Of  Fuel 

M1 

3iiOO  K 

3100  K 

3U00  K 

D 

p°  •  .12  atm 

Hi  atm 

7.5  atm 

1.0  atm 

B 

111  atm 

16.5  atm 

8.5  atra 

11  atm 

B 

16  atm 

18  atm 

10  atm 

ij  atn 

2.2 

19  atm 

22  atm 

12.5  atm 

111  .5  atm 

.  — 
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SECTION  III 

ISOTHERMAL  DESIGN  OF  A  DCW  GENERATOR 

After  selecting  the  pressure  level  as  discussed  In  the  previous  sec¬ 
tion,  the  design  philosophy  calls  for  determining  the  Isothermal  duct 
shape,  l.e.,  the  Isothermal  core  area  variation  In  the  flow  direction. 

The  result  will  serve  as  the  base  shape  upon  which  the  wall  effects  will 
be  added  In  the  computer  simulation  that  follows  the  preliminary  design. 
To  this  end,  let  us  first  collect  all  the  governing  equations  (referring 
to  Figure  2  which  describe  a  general  core  flow: 


p  uA  =  constant  (=  p^Aj) 

(1) 

P  uu'  =  -p1  +  jyB 

(2) 

P  u{CpT'  +  (u2/2) 1 }  =  jxEx  +  jyEy 

(3) 

P  =  P/RT 

(4) 

E  =  aE 
x  y 

(5) 

K  =  Ey/uB 

(6) 

Jx  =  {<x/(l  +  0Z)}(uB){aK  -  8(K  -  1)} 

(7) 

jy  «  {o/(l  +  02)}(uB){agK  +  (K  -  1)} 

(8) 

where  prime  denotes  differentiation  with  respect  to  x;  subscripts  x  and 
y  Indicate  specific  components;  p,  A,  j,  E,  R,  cp,  and  K  are  respectively 
the  plasma  density,  cross-sectional  area  of  flow  passage,  current  density, 
electric  field,  gas  constant  of  the  plasma,  specific  heat  capacity  at 
constant  pressure  of  plasma,  and  the  loading  factor;  all  quantities  are 
local.  In  Equation  1,  p  and  u1  are  known  from  the  procedure  presented 
In  the  previous  section;  A^  can  also  be  decided  roughly  on  the  power 
desired  as  explained  at  the  end  of  the  previous  section. 

In  the  above.  Equations  2  through  8  can  be  easily  combined  Into  the 
following  key  general  equation: 

(p/RT)c  T'  «  p'  +  {o/(l  +  62)}(uB2){a2K2  +  (K  -  1)2}  (9) 


Figure  2.  A  Duct  with  Diagonal  Conducting  Wall 
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Before  we  narrow  down  our  scope  to  the  Isothermal  case,  let  us  point 
out,  In  passing,  one  Interesting  fact  about  the  general  Equation  9:  for 
the  case  where  a  =  -8  and  K  =  0.5,  It  reduces  to  the  form 

(p/RT)cpT‘  =  p*  +  <j(uB2/4)  (10) 

which,  being  exactly  the  case  for  8  =  0  (and  K  =  0.5),  has  been  exten¬ 
sively  studied  for  all  values  of  K0  In  the  literature  (see  Reference  6 
for  a  unifying  approach).  The  case  with  a  =  -8  and  K  =  0.5  Is  practic¬ 
ally  of  great  Interest  since  a  *  -8  maximizes  the  power  density  with 
respect  to  a  ,  while  K  «  0.5  Is  not  far  from  being  optimal  unless  8 
deviates  very  much  from  1.  For  such  a  case,  all  the  available  solutions 
In  the  absence  of  the  Hall  effect  can  be  directly  employed  by  substi¬ 
tuting  0.5  for  K. 

Going  back  now  to  the  Isothermal  case  exclusively,  we  have  from 
Equation  9 


P'  *  -{o/(l  +  82)}(uB2)Ca2K2  +  (K  -  l)2)  (11) 

But,  from  Equation  2  and  8, 

{o/(l  +  82)HuB2)  =»  { p  uu 1  +  p')/{a8K  +  (K  -  1)}  (12) 


Thus,  combining  Equations  11  and  12,  we  have 


p'  =  -{a2K2  +  (K  -  1 )2}{ p  uu1  +  p'}/{a8K  +  (K  -  1  )}(13) 
or 

K{K(a2  +1)  +  (a8  -  1)>P'  -  -(a2K2  +  (K  -  1 )2}(p/RT)(uu' ) 
Integrating,  we  obtain  finally  the  relation  for  constant  a,  8,  and  K: 


p/pi  =  exp  {a(u12  -  u2)} 

where 

a  .  a2*2  +  (K  -  l)2 

(2RT)K{K(a2  +  1)  ♦  (a&  -  1)} 


(14) 

(15) 
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For  a  variable  g 


*  b/p,  but  with  a  =  -g  and  constant  K,  we  have 


(u^  -  u2)/2RT  *  In 


K  -  1 

P_\K  /pf(K  -  l)2-!-  bK 
P1  /  \p2(K  -  1)2+  bK 


2K  -  1 
2K(l  -  K) 


After  establishing  such  a  p  vs.  u  relationship,  we  can  go  back  to 
Equation  2  and  Integrate  It  with  respect  to  x.  The  result,  with  the 
help  of  Equation  13,  Is  as  follows: 


x 


(1  +  82)dp 

{a2K2  +  (K  -  l)2}  (auB2) 


07) 


where  the  p  vs.  u  relationship  Is  to  be  Incurred  In  the  Integration. 

Note  also  that  Equation  17  Is  completely  general,  with  varying  a,  g,  etc. 


In  any  given  case,  p(x)  obviously  will  become  known  at  this  step; 
u(x)  will  then  emerge  from  the  p  vs.  u  relationship.  Next,  since  p*p, 
A(x)  can  be  calculated  via  Equation  1.  To  determine  the  total  length 
of  the  generator,  one  sums  up  the  power  output  along  the  duct. 

Although  It  Is  not  likely  that  the  exit  pressure  would  turn  out 
lower  than  0.2  atm,  this  value  must  be  checked,  since  exit  pressure 
lower  than  0.2  atm  will  burden  the  diffuser  heavily.  If  the  magnetic 
field  Is  high  enough  that  the  foregoing  procedures  yield  a  pressure 
level  that  cause  diffuser  loading  problems,  then  engineering  compromises 
must  be  made  between  diffuser  design  and  magnetic  field  strength  In 
order  to  end  with  a  feasible  design. 


Based  on  such  a  preliminary  core  design,  end  effects,  wall  losses, 
as  well  as  real-gas  effects  can  be  added  In  computer  models  to  yield 
the  adjusted  duct  shape.  It  Is  not  known  whether  the  above  "almost- 
Isothermal"  procedure  Is  actually  being  used  by  the  designers  In  the 


13 


AFWAL-TR-80-2077 


field.  But  in  at  least  two  real  cases,  both  for  supersonic  and  for  high 
power  density,  the  temperature  deviation  Is  below  2.5%:  It  decreases  by 
about  2.3%  In  the  AVCO  400  kW  generator  (Reference  7),  and  apparently 
Increases  by  1.3%  In  the  Soviet  PAMIR  I  generator  (private  communication). 
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SECTION  IV 
SCALING  LAWS 

The  word  "scaling"  as  used  here  Is  restricted  to  mean  dynamic  or 
geometric  similitudes  among  a  series  of  MHD  generators.  The  aim  of 
scaling  Is  to  ensure  that,  once  one  member  of  the  series  is  probed 
extensively,  either  by  actual  measurements  If  It  Is  already  operating, 
or  by  detailed  numerical  simulation  If  It  Is  not  yet  built,  the 
performance  of  any  other  member  of  the  series  can  be  predicted. 


The  fundamental  tool  behind  scaling  Is  dimensional  analysis  (or, 
alternatively,  ordering).  In  principle,  every  physical  phenomenon  can 
be  described  by  a  list  of  dimensionless  modeling  parameters.  Two 
specific  Instances  of  that  phenomenon  are  mutually  transformable  in  a 
dynamically  similar  manner  If  they  have  the  same  numerical  values  for 
all  the  dynamic  modeling  parameters.  Scaling,  or  modeling  (see  Chap¬ 
ter  4  of  Reference  3),  calls  for  keeping  the  dimensionless  parameters 
the  same  among  Individual  members  of  a  series. 

In  practice.  If  the  dimensionless  parameters  governing  a  certain 
phenomenon  are  many,  a  dynamically  similar  series  may  very  well  end  up 
containing  only  one  member.  The  conclusion  of  the  scaling  process  then 
becomes  degenerate  and  trivial— namely,  every  specific  Instance  Is  only 
similar  to  Itself.  This  would  be  the  end  of  scaling,  requiring  Indivi¬ 
dual  Instances  be  separately  probed  numerically  or  experimentally. 

In  order  to  have  non-trlvlal  scaling,  certain  dimensionless  para¬ 
meters  may  have  to  be  allowed  to  vary  from  member  to  member.  If  the 
parameters  left  open  in  this  manner  are  the  geometric  ratios  and  angles, 
we  have  a  distorted  scaling;  otherwise,  we  have  a  partial  scaling. 

For  MHD  generators,  the  scaling  will  be  severely  distorted  and  will 
require  numerical  simulation  to  compensate  for  the  drastic  deviations 
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In  scaled  design.  Since  scaling  Is  basically  a  tool  used  to  provide 
quick  guidelines  before  detailed  analysis,  a  distorted  scaling  Is  of 
no  use  to  MHO  generator  design— one  might  as  well  skip  to  numerical 
simulation  directly. 


The  rest  of  this  section,  therefore.  Is  devoted  exclusively  to  the 
partial  scaling  of  MHD  generators. 

Applying  the  ordering  process  (see  Chapter  5  of  Reference  2)  to 
Equations  1  through  8  (and  thereby  Ignoring  end  and  wall  effects),  we 
Introduce  the  following  representative  quantities: 

Density - 

Velocity - Uj 

Area™ 

x-coordlnate -  Generator  length  L  (which  Is  proportional 

to  because  of  geometric  similitude) 

Magnetic  field—  Peak  magnetic  field  B* 

Pressure——  p^ 

Temperature——  T^ 

Conductivity—— 

Current  density -  ct.UjB* 

Electrical  field— —u^B* 

Substituting  Into  these  equations  the  dimensionless  quantities  p  =  f/p^, 
etc.,  with  the  tilde  signifying  nondlmenslonallzatlon,  we  have 

p  uA  =  1 

YH12(puS')  =  -P'  +  SjyB 

p  u(T'  +  C(y  -1)/2]M12(u2/2)'}  =  [(Y-  l)/Y]S(jxEx+jyEy) 
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P  =  >/T 


K  =  Ey/uB 

3X  =  (o/O  +  32}}  (uB)  {aK  -  &(K  -  1)} 

3y  =  {a/O  +  82))  (uB)  {apK  +  (K  -  1 )} 

where  the  prime  now  denotes  differentiation  with  respect  to  (x/L).  As 
a  result  of  this  process,  a  number  of  governing  (dimensionless)  para¬ 
meters  show  up  naturally  and  unambiguously;  namely, 

S  «  o^u^L/pj 

Mi  -  Uj/ZiRtj' 

B  =  B(x/L)  *  B(x) 

0  *  o(T,  p) 

8  =  (B(T,  P))  B 

a,  K,  y 

A  series  of  geometrically  similar  generators  must  have  the  same 
numerical  values  or  variations  (for  B,  p,  and  a)  In  order  to  be  dynami¬ 
cally  similar,  as  far  as  the  core  is  concerned.  Out  of  this  list,  y  stays 
around  1.1  for  all  products  of  combustion;  we  can  therefore  omit  It 
from  further  consideration. 


From  the  literature,  we  find  that 


and 


B  *  BvT/p 

a  «  fm/pn 
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where  m  and  n  are  universal  constants  depending  on  the  temperature  range 
for  all  plasmas,  but  where  the  proportionality  constants  differ  for 
different  plasmas  .  That  Is  to  say, 

£  =  b(b*/b1)/I7/p1 

and 

5  -  i"/Pn 

So,  cr(T,  p)  retains  the  same  form,  only  the  consideration  of  Is 
needed.  Similarly,  although  B(T,  p)  retains  the  same  form  (assuming 
fixed  B(x)),  B.J  must  be  considered. 

In  practice,  a  Is  made  close  to  -8;  so,  there  Is  no  necessity  to 
Include  Its  variation.  The  loading  factor  K  is  either  kept  constant 
for  a  series  of  generators,  or  It  Is  to  vary  only  slightly;  In  addition. 
Its  Influence  on  the  power  density  can  be  estimated  using  a  slug-flow 
model.  Thus,  finally  the  list  of  governing  parameters  is  reduced  to 
S,  Mj ,  B(x),  and  8^  (with  embedded  In  S). 

These  governing  parameters  are  also  obtained,  for  example,  by 
Garrison,  Brogan,  Nolan,  et  al .  (Reference  8),  presumably  through 
standard  dimensional  analysis.  The  parameter  S  also  frequently  appears 
In  the  literature  by  way  of  dimensional  analysis;  but,  usually.  It  Is 
the  only  one  mentioned. 

Out  of  these  four  remaining  quantities,  B(x)  Is  not  likely  to  be 
fixed  for  a  series  of  generators;  we  will  have  to  start  our  partial 
scaling  by  Ignoring  the  requirement  of  a  fixed  B(x).  The  Influence 
of  the  field  strength  Is  then  Invested  only  with  the  peak  B*  which 
appears  In  S  and  In  8^  through  .  Such  a  partial  scaling  has  been 
carried  out  and  applied  by  Garrison,  Brogan,  Nolan,  et  al.  (Reference 
8).  Although  It  Is  somewhat  fruitful  and  useful,  such  scaling  yields 
rather  restrictive  modeling  laws. 
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Noting  that  the  g-level  In  a  well -designed  generator  Is  around  1, 
and  that  the  Influence  of  the  e-level  can  be  assessed  separately  (e.g.. 

In  the  manner  discussed  In  Section  II),  this  report  will  go  one  step 
further  and  recommend  a  partial  scaling  based  on  S  and  only.  By 
relaxing  the  requirement  that  e^  be  kept  constant  from  case  to  case 
In  a  series  of  generators,  more  productive  scaling  laws  are  possible; 
and  the  scope  of  application  Is  suddenly  widened.  For  Instance,  taking 

o  =  cT10/p0-5  (18) 

where  c  Is  a  coefficient  the  value  of  which  Is  available  for  specific 
plasmas,  we  see  that  S  =  constant  for  different  cases  In  a  series  Implies 
that1 

L  -  p1,5B-2T-l 0.5/c  (19) 

where  M  *  constant,  l.e. , 

_0.5  .... 

u  «  T  (20) 

(noting  that  the  molecular  masses  for  all  products  of  combustion  stay 
roughly  around  35)  Is  also  enforced.  Now,  the  object  In  building  a  MHD 
generator  Is  to  realize  an  electrical  power  output  W  which  can  be 
calculated  by  Integrating  the  solutions  of  Equations  1  through  8  with 
respect  to  x.  By  examining  the  dimensionless  forms  of  these  equations, 
we  conclude  that  W  In  the  form  of  a  dimensionless  parameter  must  come 
out  a  function  of  (under  the  present  partial  scaling)  S  and  M,  without 
the  necessity  of  actually  solving  the  equations.  To  be  more  specific, 
we  conclude  that 

W/(ou2B2L3)  =  F(S,  M) 

which  yields,  since  S  and  M  are  both  kept  fixed, 

W/(ou2B2L3)  =  constant 

l  From  this  point  on,  the  subscript  1  and  superscript  *  will  be  omitted 
for  ease  of  writing. 
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W  *  ou2B2L3 

«  cT11B2L3/p0-5  (21) 

or,  In  terms  of  the  power  density, 

w  cc  cT^V/p0,5  (22) 


Next,  combining  Equations  19  through  22  In  various  ways,  we  have 
the  following  additional  scaling  laws: 


L  <*  W 

p  cc  c 

w  «  B 


0.375/(c0.25t2.75b0.5 

0.5t6.6w0.25b 

1.5c0.75t7.7 


) 


(23) 

(24) 

(25) 


1  1 25  3 

v:here  Equation  25  is  approximate  in  the  sense  that  W  /L  ,  Instead 
of  W/L3,  Is  regarded  as  being  proportional  to  w.  Furthermore,  Equation 
23  squared,  multiplied  by  Equation  24,  yields 

W-pL2/!-1*1  (26) 


In  applying  these  scaling  laws,  one  must  always  bear  In  mind  four 
things:  (1)  A(x)  Is  fixed;  (2)  Equations  18  through  20  must  be  satis¬ 
fied;  (3)  there  are  possible  errors  In  Ignoring  B(x),  a,  8,  and  K;  and 
(4)  deviations  must  be  anticipated  from  end  and  wall  losses.  To  guard 
against  possible  misapplications,  let  us  quote  here  one  counter-example: 
In  applying  Equation  26  to  the  same  generator,  one  sees  that  W  Is  pro¬ 
portional  to  p;  but  this  Is  singularly  uninteresting,  since  Equation  19 
then  dictates  that  p  be  fixed  (for  fixed  B,  T,  and  c).  A  correct 
Interpretation  of  Equation  26  would  be,  for  example,  thus:  for  a  given 
generator,  equipped  with  a  different  magnet,  the  pressure  level  must  be 
adjusted  according  to  Equation  19  In  order  to  operate  In  a  dynamically 
similar  manner  as  when  the  old  magnet  Is  used;  then  (and  only  then). 
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W  will  change  In  direct  proportion  to  p.  (If  nothing  but  p  Is  changed, 
the  new  operation  will  not  be  dynamically  similar;  and  W  will  not  be 
proportional  to  p.)  We  would  also  like  to  emphasize  that  It  Is  the 
dimensionless  quotient  w/ou^B^  that  remains  fixed  from  member  to  member 
In  a  series  of  dynamically  similar  generators;  the  power  density  w 
definitely  will  change  from  case  to  case.  It  Is  especially  Important 
to  bear  In  mind  that.  If  a  given  generator  Is  designed  for  maximum  power 
density,  and  If  a  scaled-down  unit  Is  built  according  to  these  scaling 
laws,  the  smaller  unit  Is  going  to  yield  lower  power  density. 

Among  the  previously  quoted  scaling  laws.  Equations  19,  23,  24,  and 
26  have  been  derived  before  by  Rosa  (Chapter  7  of  Reference  4),  using 
an  Intuitive,  heuristic,  and  simplistic  argument  (the  requirement  of 
fixed  M  being  absent). 

In  the  rest  of  this  section,  we  will  apply  some  of  the  scaling  laws 
to  scale  up  or  down  some  existing  designs  known  for  their  high  power 
densities.  We  will  use  subscripts  1  and  2,  respectively,  to  denote  the 
base  design  and  the  scaled  unit.  Temperature  and  pressure  In  the 
combustion  chamber  are  used  In  the  calculations;  for  fixed  Mach  number, 
they  are  proportional  to  the  temperature  and  pressure  at  the  generator 
Inlet.  The  two  specific  formulas  used  are  Equations  19  and  22  which 
are  rewritten  as 


L2/Ll 

-  (c1/c2)(p2/p1 )] •5(B1/B2)2(T1/T2)10*5 

(27) 

=  (c^KT^)11^^)0-5 

(28) 

The  results  of  the  scaling  are  summarized  In  Table  2.  A  detailed  descrip¬ 
tion  of  the  base  designs  Is  given  In  the  following: 

Base  Design  #1— AFAPL  KIVA-1  (Reference  9) 

Toluene  +  02,  seeded  with  Cs 
p°  -  10  atm,  T°  -  3100  K 


21 


AFWAL-TR-80-2077 


Mi  *  2,  B*  »  2.3  T 
W  *  200  kW,  w  *  40  MW/m3 
L  *  0.7  m 

Inlet  dimension  24.9  mm  x  99.8  mm 
Exit  dimension  72.6  mm  x  114.3  mm 

Base  Design  #2— Soviet  PAMIR-1  (Private  communication) 

Solid  fuel  (c  -  2.22xl0“33  (mho/m)(atm)°*5(K)-10) 
p°  *  45  atm,  T#  *  3559  K 
Mf  »  2.14,  B*MT 
W  3  15  MW,  w  ■  500  MW/m3 
L  *  1  m 

Inlet  dimension  160  mm  x  140  mm 
Exit  dimension  160  mm  x  220  mm 

Base  Design  #3— Maxwell  30  MW  Unit  (Reference  10) 

JP-4  +  O2,  seeded  with  Cs 

(c  *  8.45x1 0-3^  (mho/m) (atm)0, 5  (K)-10) 
p°  =  30  atm,  T°  «  3530  K 
M1  =  2.2,  B*  -  4  T 
W  =  30  MW,  w  =  200  MW/m3 
L  3  1.3  m 

Inlet  dimension  200  mm  x  200  mm 
Exit  dimension  450  mm  x  450  mm 

Base  Design  #4— AVC0  VIKING-1  (Reference  11) 

Toluene  +  02,  seeded  with  Cs 
p°  -  15  atm,  T°  -  3100  K 
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M1  *  2.2,  B*  =  2.8  T 
W  =  2  MW,  w  =  47  MW/ra3 
L  *  1.75  m 

Inlet  dimension  50  mm  x  150  mm 
Exit  dimension  166  mm  x  249  mm 

In  the  calculation,  the  scaled  unit  always  uses  toluene  +  02, 
seeded  with  Cs,  with  c  ■  1 . 68x1 0-33  (mho/m)(atm)°‘5(K)  10. 
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TABLE  2 

DESIGN  VIA  SCALING  LAWS 
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SECTION  V 
LOSSES 


The  losses  near  the  walls  and  ends,  which  are  totally  Ignored  In 
the  foregoing  sections,  will  eventually  be  accommodated  In  the  computer 
simulation  that  must  Intercede  between  the  preliminary  design  and  the 
actual  construction.  However,  a  few  qualitative  remarks  about  wall 
effects  are  In  order. 

First  of  all,  since  MHO  generation  Is  a  volume  phenomenon,  while 
wall  losses  are  surface  phenomena,  the  percentage  loss  out  of  the  total 
power  must  be  proportional  to 

^Surface  of  the  generator)  a  1 
Volume  of  the  generator)  L 


Thus,  a  larger  unit  (scaling  up)  will  suffer  less  from  the  wall  losses, 
and  will  have  a  larger  efficiency.  (The  same  trend  Is  to  be  expected 
also  regarding  the  end  losses,  since  the  end  regions  will  occupy  a  smaller 
percentage  for  larger  units.) 


Secondly,  there  Is  the  Reynolds  number  as  a  measure  of  the  viscous 
effect,  which  Is  completely  Ignored  In  the  preliminary  design.  As  a  rule 
of  thumb,  the  scaling  should  not  change  the  Reynolds  number  by  a  factor 
of  more  than  3.  Obeying  this  rule,  one  Is  usually  sure  that  the  viscous 
effects  will  change  only  quantitatively  by  a  rather  slight  degree  (e.g., 
the  friction  coefficient  Is  roughly  proportional  to  the  one-fifth  power 
of  the  reciprocal  of  the  Reynolds  number;  also  the  disturbance  due  to 
wall  roughness  becomes  less  for  larger  units).  But,  a  drastic  change 
In  Reynolds  number  must  always  be  scrutinized  carefully,  for  fear  that 
some  qualitative  deviation  may  evolve;  the  boundary  layer  may  separate 
from  the  wall,  for  Instance. 
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